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Abstract. The contributi on presents an analysis of a ra te- independent non-local damage 
model, recently proposed by lMielke and Roubfcek (|2006l) . An analytical as well as numer- 
ical solution of a simple one-dimensional bifurcation problem is performed, demonstrating 
that, for the elementary localization test, the model is free of pathological features. 
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1. Introduction 

In engineering mechanics, damage is understood as a load-induced evolution of mi- 
crostructural defects, resulting in a reduced macroscopic material integrity. The phe- 
nomenological constitutive models of damage incorporate the irreversible phenomena by 
reducing the secant modulu s of elast i city d epending on the internal damage variable. Since 
the seminal contribution of iBazantI fll976l ). it has been well-understood that such descrip- 
tion within the framework of local (i.e. scale-free) continuum mechanics leads to an ill- 
posed problem, resulting in localization of damage growth into an arbitrarily small region. 
As a remedy to this pathology, a plethora of non-local rate-independent continuum the- 
ories, based on integral, explicit and implicit gradient approac hes, has been proposed to 
introduce an internal length scale into the description, see e.g. (iJirasek and Bazantl . 
Chapter 26) for a representative overview. 

Despite a significant increase in objectivity offered by the enhanced continuum theo- 
ries, the non-local damage formulations often suffer from the fact that the non-local vari- 
ables are introduced into the model in an ad-hoc fashion, thus violating basic constraints 
of thermodynamics. In addition, due to violation of the principle of local action, such 
incor isistencies are rather difficult to detect, especially in the m ulti-dirn ensional setting, 
e.g. flPeerlings et al.l . 120041 ). Fortunately, as first d emonstrated bv Jirasekl and later 

confi r med by a number of irideperident studies, e.g . ( Peerlings et al.l . 1200 ll ; Ijirasek and Rolshovenl . 
2003 : Di Luzio and Bazant . l200,4 lEngelen et al.l . l2006l ) . a simple one-dimensional study of 
the localization behavior can serve as a convenient "filter" test, allowing to pinpoint various 
inconsistencies in the constitutive model formulation. 

Motiv ation of the present work ar ose from an energetic non-local damage model pro- 
posed bv lMielke and Roubicekl (120061). which conib ines the basic features of an engineering 
damage model due to iFremond and Nedjar with recent advances m the mathemat- 

ical theory of rate-independent irreversible processes (jMielkd . l2005l ). Such a connection 
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provides a powerful modeling and analysis framewor k, allowing at the same time a rigor- 
ous mathematical treatment of the cor nplete damag e ( Bouchitte et aP . 2007 ). theoretically 
supported numerical implementation (IMielke et al.l . 120071 ) as well as a thermodynamically 
consistent variational formulation of the non-local damage evolution problem. In this con- 
tribution, a localization study of the model will be performed to examine its qualitative 
behavior from the engineering viewpoint. Numerical as well as analytical results are pro- 
vided to illustrate its basic features with emphasis on the proper representation of the 
damage process from its initiation until complete failure. 



2. Global energetic formulation 

Let us consider a prismatic bar of length L, subjected to displacement-controlled uniaxial 
tensile loading. In the sequel, the bar will be represented by the interval Q = (— L/2; L/2), 
with boundary F (consisting of two points) subjected to the Dirichlet loading UB(t) (see 
Figure [1]), where t denotes the (pseudo-) time taken from interval (0;T). 
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Figure 1. Bar under uniaxial displacement-controlled tension. 



Following the standard thermodynamic approach to constitutive modeling, cf. (iJirasek and Bazant 



20021 . Chapter 25), a state of the system is described using an admissible displacement and 
damages fields u and cD, respectively. Formally, we write 

ueWi'\Q),u{x)\r = ur,{t)] (1) 

(2) 



u G K{t) = 



where K denotes the set of kinematically admissible displacements, Z stands for the set 

(2) 

of admissible damage values, an d (^) is the Sobolev space of functions with square- 
integrable weak derivatives, e.g. f Rektorys . 19821 ). often denoted as H^{fl). 

Within the adopted global energetic framework (IMielkd . l2005l ). the constitutive descrip- 
tion of the damage model is provided by the globally stored energy functional 



S(t, u, u) 



1 



Q{x))E (m'(x))' dx + I / Qff iQ'{x)y dx 
and by the global dissipation distance 

Qf (cD^(a;) — u}^{x)^ dx if uj'^{x) > uj^{x) a.e. in 

I 

+00 otherwise 



(3) 



(4) 



Physically, S represents the energy reversibly stored in the system, and V is the energy 
dissipated by changing the damage field from to cD^. In Eqs. and (jlj), E [Pa] denotes 
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the Young modulus of the material, Qf [Jm~^] is the amount of energy needed to disintegrate 
a unit volume of the material, and i [m] is a characteristic length of the model. Later it 
will become clear that the "+oo" term appearing in Eq. (jlj) ensures the irreversibility of 
the damage evolution; i.e. at any point, the damage variable cannot decrease in time. 

Now, given the Dirichlet loading UD^t), energetic functionals S and V and initial data 
u{0) and a;(0), the e nergeti c solut ion of the damage problem is provided by functions u(t) 
and uj{t) satisfying flMielkd . [2005h : 



Global stability: for all t G (0; T], m G K{t) and Q e Z 

£{t, u{t),uj{t)) < £{t, u{t),iu{t)) + V{uj{t),iu) (5) 
Energy equality: for all t E (0; T] 

S{t,u{t),uj{t)) + V{uj{0),uj{t)) = S{0,u{0),uj{0)) + [ V{s,u{s),uj{s))ds (6) 

Jo 

where V denotes the power of external loading. 

Although the previous two conditions present the formal definition of the energetic so- 
lutio n, the analysis itself will be performed using the time discretization technique, see 
e.g. flRektorvd . flOSl ) for a nice exposition. To that end, we introduce a uniform partition- 



ing of the time interval = to < ti = to + At < . . . < t^ = T and inductively solve the 
minimization problem 

{u{tk), uj (tk)) e Aig min [S{tk,u,Q) + V{uj{tk-i),Q)] for A; = 1, 2, . . . , (7) 

(u,(2)eK(tfc)xZ 

Note that the previous problem is independent from At, which is consistent with the as- 
sur ned rate-independent c haracter of the damage process. The theoretical results gathered 
in (iMielke and Roubicekl . l2006l : iBouchitte et all . l2007l : iMielke et all l2007f ) show that, un- 



der reasonable data qualification, the solution of the time-discretized problem converge to 
the energetic solution as At 0. 

3. Analysis of localization behavior 

Having established the essentials of the global energetic framework, we will now proceed 
with the main goal of this contribution - the analysis of the simple uniaxial tensile test. 
Intuitively, we expect that as long as the energy density at all points of the structure stays 
below Qf, the response of the bar remains elastic and no damage evolution occurs within 
the structure. After the energy density reaches the value of at a certain time tk, damage 
is initiated and localization occurs. Therefore, the explicit description of the localization 
phenomenon will build on the incremental variational principle (^^. Note that, for the sake 
of brevity, the time instant tk and the time increment At will be omitted in the rest of this 
paper. 

Dropping the constraints on u for a moment, the optimality conditions for the mini- 
mization problem ([7]) read 



Su'{x){l — uj{x))Eu'{x) — \5uo{x)E {u{x))^ + 5uj'{x)gfC^uj'{x) + 5uo{x)gf 



dx > 
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for all admissible variations 6u vanishing on T and arbitrary 6(jJ. Assuming the involved 
fields are smooth enough, we formally obtain 



- f 5u{x)(^{l-iu{x))Eu'{x)y dx- I 5uj{x) (^\E {u'{x)f + gfl^u"{x) - gj^ 



dx + 



+ [5u{x)gfeu\x)\^ > 

Using the arbitrariness of the variations and imposing the con straints < ci; < 1 finally 
leads to the system of the Karush-Kuhn- Tucker conditions, cf. ( Jirasek and Bazantl . 120021 . 
Section 15.2.4) 



({l-uj{x))Eu{x)j = (8) 

gffuj"{x) + lE{u'{x)f < gj for uj{x) = (9) 

gffuj"{x) + lE{u\x)f = gj for < uj{x) < 1 (10) 

gfe^uj"{x) + lE{u'{x))^ > gj for uj{x) = 1 (11) 

Eq. ([8]) corresponds to the uniform stress condition 

{1 - uj{x)) Eu'{x) = const = a (12) 



while Eqs. (19]), ( jlOl) and (jlTl) provide the link between strain and damage and describe 
the elastic, damaging and fully damaged part of the structure, respectively. Note that by 
virtue of Eq. (fT2l) . the stress a remains to be well-defined even in the fully damaged region, 
where the physical meaning of displacements u and strains u' is rather questionable. This 
opei i the way to the rigorous mathematical tr eatment of the complete damage problem, 
see ( Bouchitte et al. . 2007 : Mielke et al.l . 2007 ) for further discussion. 

Concentrating on the damaging zone only, we divide Eq. ffTUl) hj gj and, employing 
identity ffT^ . rewrite it in the form of a non-linear ordinary differential equation 

where 

parameterizes the softening process and ranges from /i = 1 at the onset of localization to 
/i = at the stress-free state of complete failure. It is also useful to rewrite Eq. ( |T3|) in the 
rate form 

+ 2 ^ dj{x) = - \ (15) 
(1 — uj[x)y (1 — ul)[x)Y 

which is an ordinary second-order differential equation for the damage rate io with possibly 
non-constant coefficients. Note that damage irreversibility requires ci; > 0. 
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3.1. Bifurcation from uniform state at damage initiation. Observe that at the onset 
of damage, uj{x) = and /i = 1. Therefore, Eq. ( |T5l) simphfies to 



tdj"{x) + 2dj{x) 



leading to the general solution 



ci;(x) = — |/i + Ci cos 



V2{x - 



(16) 



(17) 



with Ci and Xs denoting the integration constants, the latter specifying the position of the 
center of the localized profile. Assuming that this zone corresponds to an interval {xa, Xb) 
located sufficiently far from the boundary, the damage rate has to verify four independent 
conditions 



yielding 



Xg 



Uj{Xa) = dj{Xb) = Uj'{Xa) = dj'{Xb) = 



(Xa + Xb), Xb - Xa = Li = V27l£, Ci 



(18) 



(19) 



where Li is the width of the localized zone, cf. Figure The rate of the damage variable 
at the onset of damage is therefore given by 



UJ[XI 



— |/i 1 + cos 



y/2{x - Xs) 



7ti 7li\ 



(20) 



The previous expression is consistent with the constraint uj > since the stress rate a is 
negative in the softening regime, which implies /i < 0. 








1 ; 






— ^ ^ 



(a) -Li/2 Li/2 x ~ x, (b) -Lo/2 Lo/2 x ~ x. 

Figure 2. Sketch of bifurcation shapes; (a) initial (/i = 1), (b) final (/i = 0). 
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3.2. Complete damage. To gain insight into the behavior at complete failure, observe 
that Eq. (fT3l) . formally written for = 0, leads to a differential equation 

eu"{x) = 1 (21) 

which admits solutions in the form 



ui{x) 



1 ^jC iC <s 



2P 



(22) 



Analysis analogous to the treatment of damage initiation leads to the integration con- 
stants 

Lo = 2v^£, C2 = ^, C3 = l (23) 

resulting in the bifurcation shape shown in Figure [2b. Clearly, such a solution is not 
compatible with the irreversibility constraint, since Lq < Li. Nevertheless, it reveals the 
tendency of the damaging zone to shrink as fi decreases from 1 to 0. 

3.3. Numerical solution. The previous analysis demonstrates that the irreversibility 
condition plays an important role in the damage evolution. Therefore, the damage pro- 
file evolution needs to be determined by integrating the rate form Eq. ( fT5|) numerically, 
supplemented with irreversibility condition u > 0. As the goal of the current paper is to il- 
lustrate the basic feature of the model rather then to develop an efficient numerical scheme, 
the resulting problem is resolved by the implicit Euler ( or Rothe-Rektorys) scherne with 
the damage rate u determined using the bvp4c solver (IKierzenka and Shampind . l200ll ) 
implemented in Matlab®. The corresponding damage profile evolution, appearing in 
Figure [3^, corresponds to = 200 time steps and spatial discretization of the interval 
{xs, Xs + Li/2) into 200 elements. Moreover, the associated load-displacement diagram ap- 
pears in Figure to demonstrate that the model can properly describe complete failure. 
Both results confirm a physically admissible response of the studied model, free of locking 
effects. 
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Figure 3. Localization study; (a) evolution of damage profile, (b) load- 
displacement diagram. 
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